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Abstract 

We investigate a scenario that the top quark is the only window to the dark matter particle. We 
use the effective Lagrangian approach to write down the interaction between the top quark and 
the dark matter particle. Requiring the dark matter satisfying the relic density we obtain the size 
of the effective interaction. We show that the scenario can be made consistent with the direct and 
indirect detection experiments by adjusting the size of the effective coupling. Finally, we calculate 
the production cross section for ti + XX a t the Large Hadron Collider (LHC), which will give rise 
to an interesting signature of a top-pair plus large missing energy. 
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I. INTRODUCTION 



The presence of cold dark matter (CDM) in our Universe is now well established by a 
number of observational experiments, especially the very precise measurement of the cosmic 
microwave background radiation in the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment [lj. The measured value of the CDM relic density is 

ttcDMh 2 = 0.1099 ±0.0062 , 

where h is the Hubble constant in units of 100 km/Mpc/s. Though the gravitation nature 
of the dark matter is established, we know almost nothing about the particle nature, except 
that it is, to a high extent, electrically neutral. 

One of the most appealing and natural CDM particle candidates is weakly-interacting 
massive particle (WIMP). It is a coincidence that if the dark matter is produced thermally 
in the early Universe, the required annihilation cross section is right at the order of weak 
interaction. The relation between the relic density and the thermal annihilation cross section 
can be given by the following simple formula j^J 

Olpb 

x ~ (av) ' 1 j 

where (av) is the annihilation rate of the dark matter around the time of freeze-out. Given 
the measured f2cDM^ 2 the annihilation rate is about 1 pb or 10~ 26 cm 3 s~ 1 . This is exactly 
the size of the cross sections that one expects from a weak interaction process and that would 
give a large to moderate production rate at the LHC. In general, production of dark matter 
at the LHC would give rise to a large missing energy. Thus, the anticipated signature in the 
final state is high-p T jets or leptons plus a large missing energy. Note that there could be 
non-thermal sources for the dark matter, such as decay from exotic relics like moduli fields, 
cosmic strings, etc. In such cases, the annihilation rate in Eq. (P) can be larger than the 
value quoted above. 

The most studied dark matter candidate is perhaps the neutralino of the supersymmetric 
models with i?-parity conservation. In this work, we study a different scenario. The dark 
matter is in a hidden sector and the only standard model (SM) particle that it interacts 
with is the top quark. The top quark, having a mass so close to the electroweak symmetry 
breaking scale, makes itself so unique among the fermions. It is perhaps one of the best 



windows to probe the electroweak symmetry breaking. The dark matter, if it is a WIMP, is 
also closely related to electroweak symmetry breaking. The logic is that since both the top 
quark and the WIMP are closely related to electroweak symmetry breaking, we argue that 
the top quark may be the only window to probe the dark matter. This is our motivation, [j] 
We use an effective Lagrangian approach to parameterize the interactions between the top 
quark and the dark matter particle, without specifying the detailed communication between 
the top quark and the hidden sector. One simple example would be a hidden-sector gauge 
boson that can couple to the top quark. If it is heavy enough we can shrink the propagator 
to a 4-fermion vertex. One form of the 4-fermion interaction is (£7^) (x7 M x) f° r a vector- 
type interaction or (ij^t) (x7 M 7 5 x) for an axial-vector-type interaction. We can estimate 
the size of the new interaction based on the fact that it is the only interaction that can 
thermalize the dark matter particle in the early Universe. The most interesting implication 
of the scenario is the collider signature. The final state consists of a top-quark pair and a 
pair of dark matter particles, giving rise to a top-quark pair plus a large missing energy. On 
the other hand, we anticipate that the spin-independent and spin-dependent cross sections 
in direct detection would be consistent with the current limit. This is easy to understand 
because the top content inside the nucleon is so small that it hardly contributes to the 
DM-nucleon scattering. We will explicitly show that. In addition, the annihilation of the 
dark matter in the Galactic Halo would give rise to positrons and antiprotons that can be 
observed by antimatter search experiments, e.g., PAMELA and AMSII. We use the present 
data on positron and antiproton spectra from PAMELA to constrain the size of the effective 
interactions. 

The organization of the paper is as follows. In the next section, we describe the interaction 
between the top quark and the dark matter particle, and estimate the size of the interaction 
based on the relic density. In Sec. Ill, we calculate the spin-independent and spin-dependent 
cross sections for direct detection. In Sec. IV, we calculate the positron and antiproton 
spectra due to the DM annihilation in Galactic halo. In Sec. V, we discuss the collider 
signature. Finally, we conclude in Sec. VI. 

nn 

There are a few recent works |4J-|6| that assumed some form of effective interactions be- 
tween the dark matter and light quarks and studied the corresponding collider phenomenol- 



ogy. Fan et al. 



7| also wrote down the effective nonrelativistic interactions between the 



1 A possibility of realizing such a scenario can be found in Ref. 
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dark matter and nuclei. 



II. EFFECTIVE INTERACTIONS AND RELIC DENSITY 

Our simple model consists of the SM and a hidden sector, in which there is a pair of 
Dirac/Majorana fermions and a gauge boson. For some reasons this gauge boson couples 
this hidden fermion only to the top quark on the SM side. If the mass of this gauge boson 
is heavy enough, we can integrate it out. More generally, below the heavy mass scale A the 
interaction between the top quark and the dark matter particle, denoted by \i is given by 

£ = ft (xTx) (in) , (2) 

where T = 7^ for a vector gauge boson, T = 7^7 5 for an axial-vector gauge boson, T = 1 (7 s ) 
for a scalar (pseudoscalar) boson interaction, and V = <r MV (7 5 ) with <j^ v = i{ r y^ r y u — 7 !/ 7 At )/2 
for a tensor (axial-tensor) interaction, and g x is an effective coupling constant. For Majorana 
fermions the T = 7^ or a^ v type interaction is identically zero, and so for vector or tensor 
type interaction the fermion x lYi Eq-tEl) must be Dirac. Explicitly, we assume the dark 
matter candidate to be Dirac, but the results are also applicable to Majorana dark matter. 
With this interaction we can calculate the thermal averaged cross sections and thus the 
relic density, the direct and indirect detection rates, and also the production cross section 
of pp — > ti + XX a t t ne LHC. 

We start with a vector gauge boson type interaction: r = 7^. The differential cross 
section for xipi) xipz) ~^ t(ki) i(k 2 ), with the 4-momenta listed in the parentheses, is 

l4iM [ "" + t ™ + 2sK + m?)1 (3) 

where Nc = 3 for the top quark color, /3 t , x = (1 — 4m^ x /s) 1//2 , t m = t — m x — vr? t = 
— s(l — f3 t f3 x z)/2, u m = u — m x — mf = —s(l + (3 t (3 x z)/2, s = (pi +P2) 2 is the square of the 
center-of-mass energy, t = (p\ — &1) 2 , u = (pi — /C2) 2 , and z = cosO with the scattering 
angle. The quantity av, where v ~ 2(3 X in the non-relativistic limit, can be obtained by 
integrating over the variable z in Eq. ([3]). Instead of solving the Boltzmann equation, we 
can naively estimate the size of the interaction by the following equation 

syf (4) 

x (av) 




FIG. 1. The calculated av versus for the effective interaction ^ (x^x) (*ri) of various Dirac 
structures T with A = 1 TeV, m x = 200 GeV, and v « 0.3. 

With the most recent WMAP result on dark matter density n CDM /i 2 = 0.1099 ± 0.0062 we 
obtain the size of av 

(av) ~ 0.91 pb . (5) 



We show in Fig. [T] av versus the coefficient g 2 x for a dark matter mass of 200 GeV. The 
result shown is for v ~ 0.3 to approximate the velocity of the dark matter particle at around 
the freeze-out time. We can repeat the calculation with T = <r^(7 5 ), 7 At 7 5 , 7 5 , 1 for tensor 
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(axial-tensor), axial- vector, pseudoscalar, and scalar type interactions. The results are 



da g x N c /3 t 



dz A 4 47rs p % [ 2 (4 + O + 2 *K + O + §™K - ^ ( 6 ) 

da g x N c f3 t r 2 2 2 2 2 21 

dl = A 1 16^ F ^ m Mm ~ (m * ™ x) 16m * m J ( ? ) 



da = &Nc_ fk (8) 
dz A 4 32vr (3 x 1 J 

rfcr gi N, 



dz A 4 32vr 



sPxft (9) 



for T = a fiu ('-f 5 ), 7 M 7 5 , 7 5 , 1, respectively. We note that the axial-tensor case has the 
expression as in the tensor one given by Eq.([n]). The results are shown in Fig. [Has well. We 
can see that the tensor-type interaction gives the largest cross section, followed by vector, 
pseudoscalar, and axial-vector. These four types of interactions require g x falling into the 
range of 0.2 — 0.6 which is about the size of weak-scale interaction. On the other hand, 
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the scalar-type interaction always gives a very small annihilation cross section for a similar 
range of g x , which is in danger of over-closing the Universe. 

In Fig. [2j we show the contour of the cross section for the various types of interactions 
as a function of g x and m x as allowed by the WMAP result. 



III. DIRECT DETECTION 



a. 



Recently, the CDMSII finalized their search in Ref. [8j. When they opened the black 
box in their blind analysis, they found two candidate events, which are consistent with 
background fluctuation at a probability level of about 23%. Nevertheless, the signal is not 
conclusive. The CDMS then improves upon the upper limit on the spin-independent (SF 
cross section af N to 3.8 x 10" 44 cm 2 for m x » 70 GeV. The XENON100 Collaboration [9j 
also recently announced their newest result. Although XENON100 has the best sensitivity 
in the lower mass range, the CDMSII limit is currently still the best in the world for dark 
matter mass larger than about 100 GeV. We will adopt a limit of order 4— 10 x 10 -44 cm 2 for 
dark matter mass of 200 — 500 GeV. In the following, we will check if the spin-independent 
cross section generated by the 4-fermion interactions is consistent with the new limit. 

Spin-independent cross sections can arise from the scalar-type and vector-type interac- 
tions between the DM and quarks. If the effective interactions between the dark matter 
particle and the quarks are given by 

£= Yl i^XXQQ + ^xYXQltiQ} , (10) 

q=u,d,s,c,b,t 

then the spin-independent cross section between \ an d each of the nucleon (taking the 
average between proton and neutron) is given by 



*? N = -^ iGff + ^H , (id 



vr V s 1 256 

where fi x ^ = m x m N j (m x +m N ) is the reduced mass between the dark matter particle and 
the nucleon N, and 

Gs= E (N\m\N) <* S q , (12) 

q=u,d,s,c,b,t 

where (N\qq\N) denotes the various nucleon matrix elements. The expression for of a 
whole nucleus (A, Z) is = ot l(A + Z) + a^{2A - Z), we take the average between proton 



and neutron (assume the number of protons is about the same as that of neutrons in the 
nuclei) and thus obtain the expression for a single nucleon 

b N = 3 -(a v u +a v d ) . (13) 

Nevertheless, the contributions to come from valence quarks only. Therefore, in our 
scenario the only contribution to the SI cross section comes from the top quark, thus only 
af is nonzero in the expression of , which is then given by 

Gf = (N\tt\N) (&) , (14) 



DarkSUSY 



V A 2 

where (N\tt\N) = f^ t {mM /ra t ). □ The default value of the parameters f^ t used, e.g. in 



10|, is 

fP f = 0.0595 , fa = 0.0592 



Taking the average between proton and neutron the value of is 



G^^(ll). (15) 



m t \A 2 

For m x ~ 0(100) GeV, h x n ~ m^. The spin-independent cross section is 

^^feVfiiV. d«) 



We show in Fig. [3] the spin-independent cross section versus gi. Note that the axial- 
vector interactions contributes to spin-dependent cross sections. Since the constraint from 
spin-dependent cross sections is a few orders of magnitude weaker than that from spin- 
independent cross sections, we simply focus the spin-independent one to obtain the mean- 
ingful range of g x and A. We found that the limit on spin- independent cross section of 
the order of 10~ 44 cm 2 allows gi as large as 30 for A = 1 TeV. Note that for a strongly 
coupled theory, one can have g 2 = (47r) 2 . Such a large g x is allowed by spin-independent 
cross section constraint as well as by the WMAP relic density constraint. □ However, one 
must be cautious that for such a large effective coupling constant, perturbative calculation 
becomes less reliable. 

We next turn to the indirect detection of the dark matter, which then gives the strongest 
constraint on the present scenario. 

2 Equivalently, the top content inside the nucleon can be replaced by the gluon content with f^ q replaced 

by ^ffxg 21 ■ Numerically, they are very close to each other. 

3 When the annihilation cross section is larger than that required by thermal production, the resulting 

relic density from thermal production is just too low. However, there could be some other non-thermal 
sources, such as decay from heavier fields. 
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x 



FIG. 3. Spin-independent cross sections for the vector type interaction versus g^. 



IV. INDIRECT DETECTION 



Another important method to detect the dark matter is by measuring its annihilation 
products in Galactic halo. Current experiments can detect the positron, antiproton, gamma 
ray, and deuterium from dark matter annihilation. The Milky Way Halo may contain clumps 
of dark matter, from where the annihilation of dark matter particles may give rise to large 
enough signals, such as positron and antiproton, that can be identified by a number of 
antimatter search experiments. The most recent ones come from PAMELA 111. Il2l|. which 
showed a spectacular rise in the positron spectrum but an expected spectrum for antiproton. 
It may be due to nearby pulsars or dark matter annihilation or decays. If it is really due 
to dark matter annihilation, the dark matter would have very strange properties, because 
it only gives positrons in the final products but not antiproton. Here we adopt a conserva- 
tive approach. We use the observed antiproton and positron spectra as constraints on the 
annihilation products in XX annihilation. 

We first consider the positron coming from the process 



XX -> tt ->• {bW + )(bW-) -> (be + v e ) + X 



(17) 
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in which the most energetic e + comes from the W + decay. There are also positrons coming 
off in the subsequent decays of b,b, r + , or fi + , but these positrons are in general softer 
than those coming directly from the W + decay. For a first order estimate of the size of the 
coupling g^, in Eq. §2§ we only include the positrons coming directly from the W + decay. 

The expression for annihilation has already been given in Eq. (J3]), but now with a present 
time velocity v ~ 10 -3 . The positron flux observed at the Earth is given by 

$ e+ (E) = ^fe + (E), (18) 

with v e + close to the velocity of light c. The function f e + (E) satisfies the diffusion equation 
of 

Bf 8 

-± - K(E)V 2 f - — (b(E)f) = Q , (19) 

where the diffusion coefficient is K(E) = K (E/GeV) s and the energy loss coefficient is 
b(E) = E 2 /(GeV x te) with te = 10 16 sec. The source term Q due to the annihilation is 

\ 2 

PCDM \ V^/ \ dN e 



where rj = 1/2 (1/4) for (non-)identical DM particle in the initial state. The summation 
is over all possible channels that can produce positrons in the final state, and dN e + / dE e + 
denotes the spectrum of the positron energy per annihilation in that particular channel. In 
our analysis, we employ the vector-type interaction for Dirac fermions and thus the source 
term is given by 

n 1 ( PCDM V , v dN e + 

where the normalization of N P + is 



/ 



^^dx = B(t -)• bW + -)• be + u e ) . (22) 

We then put the source term into GALPROP [13!] to solve the diffusion equation. In Fig. H] 
we show the predicted energy spectrum for the positron fraction for various values of g^. 
With a visual inspection the < 8 is allowed by the spectrum. 

Next we turn to the antiproton fraction as it was also measured by PAMELA. Similarly, 
the antiproton flux can be obtained by solving the diffusion equation with the corresponding 
terms and the appropriate source term for the input antiproton spectrum: 

2 

M CDM ) ^""' v dTp 
10 



Q^-v(^) EH-f, (23) 



♦ pamela data 
galprop bkgd 

— g_x A 2=2 

g_x A 2=4 

g_x A 2=6 

g_x A 2=7 




1000 



energy (GeV) 



FIG. 4. Spectrum for the positron fraction predicted for the vector type interactions for various 
g^. PAMELA data are shown. 

where 77 = 1/2 (1/4) for (non-)identical initial state, and Tp is the kinetic energy of the 
antiproton which is conventionally used instead of the total energy. We again solve the 



diffusion equation using GALPROP 



13| 



In our case, the dominant contribution comes from 

XX ->it-> (bW + )(bW~) -> (bqq')(bqq') -> p + X . (24) 
In the last step, all the b b, q, q' have probabilities fragmenting into p. We adopt a publicly 



available code (14| to calculate the fragmentation function D q ^ h (z) for any quark q into 
hadrons h, e.g., p,p, n. The fragmentation function is then convoluted with energy spectrum 
dN/ dE of the light quark to obtain the energy spectrum of the antiproton dN/ dEp. The 
source term dN/ dTp is then implemented into GALPROP to calculate the propagation from 
the halo to the Earth. We display the energy spectrum for the antiproton fraction in Fig. [51 
It is easy to see that is constrained to be 

<£<4-5. (25) 
11 
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FIG. 5. Spectrum for the antiproton fraction predicted for the vector type interactions for various 
g^. PAMELA data are shown. 

We will use this allowed range to estimate what we would expect from the LHC. 



V. COLLIDER SIGNATURE 

Collider signatures are perhaps the most interesting part of the scenario - ti pair plus 
large missing energy. We first calculate using the effective 4-fermion interaction with F = 7 M 
the production cross section for pp — >■ ti + XX- There are two contributing subprocesses for 
ti production at the LHC: 

qq->it, gg->tt, (26) 

on which we can attach one 4-fermion interaction vertex to each fermion leg including internal 
fermion line to further produce a XX pah- A typical Feynman diagram is shown in Fig. |6j 



We employ MADGRAPH 15[ to calculate the signal and background cross sections. 



The irreducible background is tt + Z — > ttvv. Before applying any cuts we calculate 

12 



FIG. 6. A contributing Feynman diagram for the subprocess gg — > tt + xX- The other two diagrams 
can be obtained by attaching the black dot to the t and i leg, respectively. 



FIG. 7. Missing transverse energy fir distributions for the signal pp — > tt + xX an d the background 
pp -)• ttZ for g\ = 3 - 30 with m x = 200 GeV. 

the signal cross section versus background cross section: 8.2 fb (for = 5) to 140 fb, in 
which we have chosen scale Q = (2m t + 2m x )/2 in the running coupling constant and the 
parton distribution functions for the signal, while Q = (2m t + mz)/2 for the background. 
We first compare the missing Ej< distribution between the dark matter signal and the tiZ 
background, shown in Fig. [71 It is clear that the signal has a harder missing Et spectrum 
than the background. This plot suggests a cut as large as 400 GeV in the missing transverse 
energy can substantially reduces the background to a level similar to the signal. The cross 
sections in fb for the signal and the background using various cuts on the missing energy 
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FIG. 8. Event numbers for the invariant mass tt distributions for the signal pp — > tt + XX an d 
the background pp — > ttZ (a) before and (b) after applying the missing transverse momentum cut 
of 400 GeV. The assumed luminosity is 100 fb" 1 , which corresponds to 240 signal events and 420 
background events after the cut. 



are shown in Table [H The background can indeed be cut down to the level as the signal 
with a missing energy cut of 400 GeV. Note that the signal cross section scales as g^. The 
significance of the signal S/yB for an integrated luminosity of 100 fb" 1 stays around 11 
with a cut of 300 — 500 GeV. Since the significance scales as a/£, with a reduced luminosity 
of 30 fb" 1 the significance is still as large as 6. 

We then compare the tt invariant mass distribution between the dark matter signal 
and the ttZ background before and after applying the missing Et cuts, shown in Fig. [HI 



In Table [TTl we also show the signa 
pseudoscalar, and scalar interactions 



1 cross sections and the significance for axial-vector, 
j in decreasing order. Note that the cross section at 
the LHC for scalar interaction is not severely suppressed, in sharp contrast to the annihilation 
cross section calculated in Sec. II. 



Note that the tensor interaction is not present in current version of MADGRAPH 



13- 
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TABLE I. Cross sections in fb for the signal pp —> ti+ XX an d the background pp — > ttZ — > ti + vv 
at the LHC. We used = 5 for illustration. The signal cross section scales as g^,. The significance 
S/y/B is calculated with an integrated luminosity of 100 (30) fb -1 . 



Pt > 


pp — > 


tt + xx 


p — > ttZ — > ttvv 


S/B 




S/VB (100 (30) fb" 1 ) 


GeV 


8.2 




140.3 


0.06 




6.9 (3.8) 


300 GeV 


3.6 




10.7 


0.34 




11.0 (6.0) 


400 GeV 


2.4 




4.2 


0.57 




11.8 (6.4) 


500 GeV 


1.5 




1.9 


0.78 




10.6 (5.9) 


TABLE II. Cross sections in fb for the signal pp — >■ tt + XX f° r vector, axial-vector, pseudoscalar, 


and scalar interactions 


at the LHC 


. We have imposed the fix > 


400 


GeV cut. The S/B and 


S/y/B are shown. The background 


is from Table I. The si; 


mificance S/- 


\f~B is calculated with an 


integrated luminosity of 100 (30) fb" 


-i 












Signal cross section (fb) 




S/B 


s/Vb 


Vector 




2.4 






0.57 


11.8 (6.4) 


Axial-vector 




1.9 






0.45 


9.3 (5.1) 


Pseudoscalar 




0.82 






0.20 


4.0 (2.2) 


Scalar 




0.55 






0.13 


2.7 (1.5) 



VI. CONCLUSIONS 

In this paper we have studied an interesting scenario where the dark matter couples 
exclusively to the top quark using an effective field theory approach, with the intuition that 
both the top quark and the dark matter may be closely related to electroweak symmetry 
breaking. We did not specify any particular connector linking the SM sector and the invisible 
dark matter sector, except that this connector sector was taken to be heavy, probably at 
the TeV scale. Integrating out the heavy connector sector may give rise to effective 4- 
fermion interactions of tensor, axial-tensor, vector, axial- vector, pseudoscalar, or scalar types 
between the dark matter and the top quark. We studied the constraints of these effective 
couplings from WMAP as well as from the direct and indirect detection of dark matter at 
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CDMSII and PAMELA, respectively. 

If we require all the dark matter in the Universe comes from the thermal equilibrium, 
the coupling m 0.3 — 0.6. However, if we just require that the dark matter does not 
overdose the Universe the gt can be much larger. Since only the top quark inside the 
nucleon contributes to direct detection cross section, the coupling can be as large as 40. 
On the other hand, the strongest constraint comes from the positron and antiproton fraction 
spectra. The PAMELA antiproton spectrum constrains the coupling to be g\ < 4 — 5. 

This model can be tested at colliders with a very distinct signature, namely, tt plus 
missing energies. The top quark and antiquark would mostly have high pj- and boosted. 
The detection of such boosted top quarks has attracted some recent studies that it can 

n 

be sufficiently distinguished from the background [16j. Our results suggested that this 
interesting scenario can be testable at the LHC. 
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